arXiv: 1501.03216vl [quant-ph] 14 Jan 2015 


Sorting photon wave packets using temporal-mode interferometry based on 
multiple-stage quantum frequency conversion 


D. V. Reddy, 1 M. G. Raymer, 1 ’* and C. J. McKinstrie 2 

1 Oregon Center for Optics, Department of Physics, 

University of Oregon, Eugene, Oregon 97403, USA 
2 Applied Communication Sciences, Red Bank, New Jersey 07701, USA 
(Dated: January 15, 2015) 

All classical and quantum technologies that encode in and retrieve information from optical fields 
rely on the ability to selectively manipulate orthogonal field modes of light. Such manipulation can 
be achieved with high selectivity for polarization modes and transverse-spatial modes. For the time- 
frequency degree of freedom, this could efficiently be achieved for a limited choice of approximately 
orthogonal modes, i.e. non-overlapping bins in time or frequency. We recently proposed a method 
that surmounts the selectivity barrier for sorting arbitrary orthogonal temporal modes [Opt. Lett. 
39, 2924 (2014)] using cascaded interferometric quantum frequency conversion in nonlinear optical 
media. We call this method temporal-mode interferometry, as it has a close resemblance to the 
well-known separated-fields atomic interferometry method introduced by Ramsey. The method has 
important implications for quantum memories, quantum dense coding, quantum teleportation, and 
quantum key distribution. Here we explore the inner workings of the method in detail, and extend it 
to multiple stages with a concurrent asymptotic convergence of temporal-mode selectivity to unity. 
We also complete our analysis of pump-chirp compensation to counter pump-induced nonlinear 
phase-modulation in four-wave mixing implementations. 


I. I. INTRODUCTION 

The state of a photon is fully characterized by its he- 
licity, and the three components of its momentum. For 
guided photons or photons in beam-like geometries, the 
same degrees of freedom may be stated as polarization, 
transverse-mode profile, and energy (or frequency) [1]. 
Polarization spans a two-dimensional state space, allow¬ 
ing definitions of two-tuple orthogonal basis sets, for ex¬ 
ample perpendicular linear polarization states. Sorting 
and the application of unitary transformations on any 
choice of bases is accomplished through combinations of 
existing tools such as polarizing beam splitters, half-wave 
plates, and quarter-wave plates. Transverse modes are 
two-dimensional functions, and may be expressed in the 
basis of orbital angular momentum, transverse-spatial 
parity, or integer indices denoting eigenfunctions of some 
guiding geometry. In contrast to polarization, transverse 
modes offer a multidimensional state space, and their se¬ 
lective manipulation has been an active subject of recent 
studies [2-9]. 

For traveling (longitudinally unconstrained) fields, en¬ 
ergy (or frequency) is a continuous-variable degree of 
freedom. This spans an infinite dimensional space, and 
high-fidelity selective sorting (and by extension, any uni¬ 
tary operation) on arbitrary states in this space has tra¬ 
ditionally been limited to specific choices of basis sets. 
If the considered states have disjoint spectra, this is ac¬ 
complished using a prism or diffraction grating. That is 
the basis of wavelength-division multiplexing in optical 
communication systems. Sorting has also been accom¬ 
plished for spread-spectrum distributions of orthogonal 
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frequency subcarriers via orthognal-frequency-division- 
multiplexing (OFDM) [10], and for partially overlapping 
temporal pulses of a specific shape and temporal offset 
via orthogonal-time-division-multiplexing (OTDM) [11]. 
But orthogonal states can also be defined as superposi¬ 
tions over a continuum of frequencies with overlapping 
spectra. By Fourier transformation, such states can be 
represented in the time domain, thereby introducing the 
idea of ‘temporal modes’-a set of orthogonal functions 
of time that form a complete basis for representing an 
arbitrary state in the energy degree of freedom [1, 12]. 

Temporal modes (TMs) are a prime contender for the 
preferred basis for mapping qubits or higher-dimensional 
qudits onto photons, because their orthogonality is im¬ 
mune to chromatic dispersion and polarization rotation. 
TMs are ideal for transmission in waveguides and opti¬ 
cal fibers with no cross-talk: since all TMs at a given 
central frequency suffer the same unitary evolution un¬ 
der linear dispersion, their orthogonality is maintained 
in long-distance quantum communication through time¬ 
stationary guided media. 

The most promising approach to the sorting of light 
into a discrete set of TMs has been nonlinear wave mixing 
with pulsed pumps [13, 14]. From a general mathematical 
standpoint analogous to that in [15], it is clear that tem¬ 
poral wave-packet multiplexing is possible only by using 
a temporally nonstationary process that acts on field am¬ 
plitudes, not intensities. This is inherently true in strong- 
pump driven nonlinear wave mixing, where the weak sig¬ 
nal fields need to be synchronous with a device clock, 
defined by the arrival time of the pump pulse. Nonzero 
relative group velocities between participating fields also 
allow the device to “sample” or “measure” the complex 
temporal amplitude of the signals via a continously shift¬ 
ing inter-pulse temporal overlap, which is essential if it 
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needs to descriminate between orthogonal modes. The 
complex envelopes of the strong pump pulses provide for 
a means of programming the device to sort user-defined 
temporal modes. Such a system can be called a complex- 
pulse-envelope gate, or pulse gate for short [14]. It has 
also been called a field-orthogonal TM sorter [16]. 

It was first pointed out in [13] that quantum frequency 
conversion (QFC) by four-wave mixing (FWM) is par¬ 
tially TM-selective when using short pump pulses. The 
same is true also for QFC by three-wave mixing (TWM) 
driven by a single pump pulse, as pointed out by Eckstein 
et al, and proposed as a general method of TM sorting 
and multiplexing for quantum information applications 
[14]. The idea is to design the process so that only a 
single temporal mode among the orthogonal set has its 
frequency converted; then a simple frequency-selective 
component (diffraction grating, prism, etc.) could easily 
separate the converted wave packet from the unconverted 
ones. Recent studies furthered the concept of using opti¬ 
cal frequency conversion (FC) by nonlinear wave mixing 
for creating TM-mode-selective devices [17-20]. The role 
of FC in quantum optical communication networks has 
been reviewed in [21]. 

Ideally, the device should frequency convert a user- 
specified temporal mode component from the input field 
into the output mode in a different frequency band with¬ 
out ‘contamination’ by other temporal mode components 
of the input field in the orthogonal subspace. But it 
was found through careful modeling that the basic single- 
stage QFC process is fundamentally limited to a TM ‘se¬ 
lectivity’ of no greater than about 0.8, strongly limiting 
its usefulness [14, 17, 19, 20, 22]. High selectivity is a nec¬ 
essary condition for scalability in quantum applications, 
and a value of 0.8 falls short of error-tolerance estimates 
in both quantum communication and computing. The 
limit exists as a fundamental feature of group-velocity 
mismatched inter-pulse dynamics in a dispersive medium, 
as will be shown. We discovered that a generalization of 
the QFC schemes (either FWM or TWM) can circumvent 
this limit and achieve essentially perfect TM sorting, with 
a selectivity approaching unity [16]. The new scheme uses 
two or more spatially separated stages of QFC, with the 
effects of the stages adding coherently. We named this 
scheme temporal-mode interferometry (TMI), as it has a 
close resemblance to the separated-fields atomic interfer¬ 
ometry method introduced by Ramsey [23]. In that case, 
an atom passes through two separated laser beams. The 
first field partially excites the atom (for example, to 50% 
probability), and the second field, depending on its phase 
relative to that of the oscillating atomic dipole resulting 
from the first field, either further excites the atom (for 
example, to 100%) or fully deexcites the atom back to 
its ground state. The present case of TMI is more com¬ 
plex, in that the traveling pulse is not localized at a point 
like a typical atom, but undergoes temporal shaping as 
its frequency is being converted. Time-domain versions 
of Ramsey interference are also known in the contexts 
of photon-echo experiments and atomic-fountain clocks 


[24, 25]. 

The dependence of conversion efficiency on the rela¬ 
tive phases of the signal and control fields is what earns 
the TMI method the name interferometry. It is useful 
for real-time control and switching of various temporal 
modes. The concept of two-stage QFC, and its analogy to 
Ramsey interferometry, was subsequently independently 
developed by Clemmen et al [26], but only in the context 
of continuous-wave fields, so in this case it lacks the TM 
selectivity that is the main thrust of the present discus¬ 
sion. 

This paper details and expands the work reported in 
[16]. In particular, we present the theory of inter-stage 
mode matching. We explore the dependence of mode se¬ 
lectivity on the ratio of interaction-time and pump tem¬ 
poral duration in TWM. We provide the theory of pump- 
chirp compensation in FWM, including a family of an¬ 
alytical expressions for pump chirps. Finally, we give 
numerical demonstrations of multistage (more than two) 
TMI and the associated asymptotic gains in mode selec¬ 
tivity. 

II. II. TEMPORAL MODES 

Temporal modes are a set of complex functions that are 
orthogonal in time (or frequency), which form a complete 
basis for representing an arbitrary state in the energy 
degree of freedom [1, 12]. The state can be expressed 
as a coherent superposition of many possible frequency 
states, 

| if>j) = J dujfij((jj)a\uj)\vac) = b^\vac), (1) 

[a(w), (a/)] = 2tt6(lu - w'), \bj,b[] = S j)k , 

where af (uj) creates a monochromatic photon with fre¬ 
quency cj. By Fourier transform, the same state can be 
expressed as a coherent quantum superposition state of 
many possible ‘creation times,’ that is 

\-ipj) = J dt'ip j (t)A\t)\vac), ( 2 ) 

where (t) creates a broadband photon at time t. The 
TM mode functions ^ j(t ) fully overlap in time, but are 
orthogonal with respect to a time integral. Likewise in 
the frequency domain. The mode basis sets {fjj (t )} (or 
{Vt/(co’)}) can be any user-defined orthonormal set. The 
integer index is in principle unbound, reflecting the in¬ 
finite dimensionality of a continuous variable degree of 
freedom. For exceptionally broadband TMs, covering es¬ 
sentially the entire visible spectrum, exact orthogonality 
is not ensured [1]. But this subtlety will have negligible 
effects on practical systems where channels are defined 
with bandwidths that are narrow compared to the car¬ 
rier frequencies. 
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TM-based multiplexing has more in common with 
code-division multiplexing, as the different overlapping 
functions serve as spread-spectrum codes that mark in¬ 
formation channels. This is in contrast with conventional 
time- or frequency-based optical multiplexing, which use 
either separated short pulses or narrow spectral windows 
to define different information channels (fig. 1). 



FIG. 1. Time- and frequency-division multiplexing being con¬ 
trasted with code-division multiplexing (CDM), with different 
codes represented by time-frequency Wigner functions for or¬ 
thogonal temporal modes. The choice of the temporal-mode 
set is not unique. 


III. III. SINGLE-STAGE QUANTUM 
FREQUENCY CONVERSION 

Here, we present the theory for QFC via three-wave 
mixing (TWM) [27-31] and four-wave mixing (FWM) 
[13, 32]. To model the process, we designate the par¬ 
ticipating frequency channels by s and r for the signal 
photons, and p (and q ) for the strong pump field(s). We 
denote square-normalized classical field envelopes in the j 
channel by Aj(z , t). Quantum mechanically, Aj(z , t) rep¬ 
resent the field annihilation operators for a wave-packet 
mode of that central frequency (cjj) and envelope shape 
[1] and obey the same equations of motion as their clas¬ 
sical counterparts [13, 19]. Then the equations of motion 
for QFC are [33, 34]: 


{d z + /3pdt)Ap — i('y/2)6p \A P \ +2|j4 g | A p , 

(3a) 

(d z + p' q d t )A q = i('y/2)5 F [2 |^ p | 2 + |A/] A q , 

(3b) 

0 dz + f3' r d t )A r = ijA p A*A s 

+ iX$F \A p 2 + \A q \ 2 A r , 

(3c) 

(d z + P' s d t )A s = i^A* p A q A r 

\A P \ + \A q \ A s , 

(3d) 


where = [d/3/du;]\ UJ=U j j are the group slownesses (in¬ 
verse group velocities), and we have neglected higher- 
order dispersion, which is valid for sufficiently narrow- 
band pulses. The self- and cross-phase modulation factor 
Sp is 0 for TWM and 1 for FWM. 7 is the product of the 
effective y( 2 ) or X® nonlinearity and the square roots of 
the pump-pulse energies. For TWM, A q (z,t) = 1. The 
presence or absence of complex-conjugation of pump am¬ 
plitudes follows the convention that uj r = (d p +uj s —uj q for 
FWM, and cj r = uj p + uj s for TWM. The solution for Eq. 
(3) can be written in terms of Green functions (£,£'). 
For j G {r, 5 } and medium length /: 


Aj(l, t) = / dt' ^2 G jk (t,t')A k (0,t'). (4) 

J-00 k=r,s 

The four kernels [17] comprise a gener¬ 

alized beam-splitter transformation, representing fre¬ 
quency conversion or non-conversion of each of the two 
possible input fields of distinct carrier frequencies [35]. 
The functional forms of the Green function kernels are 
dependent on the shapes of the p (and q ) pump pulse(s). 
To determine the natural temporal modes for the prob¬ 
lem, we numerically perform a singular-value decompo¬ 
sition (SVD), also called a Schmidt decomposition, for 
each Gij(t,t'). This yields a set of singular-value func¬ 
tions (temporal Schmidt modes) with associated singular 
values (Schmidt coefficients), which we denote by p n and 
r n , which satisfy \p n \ 2 + \r n \ 2 = 1. |p n | 2 is the frequency 
conversion efficiency (CE) of the nth mode, and can be 
interpreted as the probability of QFC in the case of single 
photons. The modes are numbered in decreasing order 
of the Schmidt coefficients; the mode with the largest 
Schmidt coefficient ( i.e . the first Schmidt mode, n = 1) 
is the target TM (expected to be selectively converted). 
We denote by </> n (£) the 5-input modes, and T n (£) are 
r-output modes. In addition, there are r-input modes 
Tpn (t) and 5-output modes <£> n (£). 

All the Green function (GF) kernels can be expressed 
in terms of these four mode sets as follows [35]: 

00 

G rr (t,t') = (5a) 

n =1 
00 

G rs {tF) = Y,Pn^n{tW n {t'), (5b) 

n= 1 
00 

G sr {t, t') = - ^ (5c) 

n= 1 
00 

G ss (t,t') = ( 5d ) 

n= 1 

As mentioned earlier, the individual kernals are not 
unitary, but taken together, the complete GF is. The 
objective is to design a QFC device that can selectively 


















4 


frequency convert or ‘drop’ the first Schmidt mode (tar¬ 
get TM) with unit efficiency, whilst allowing 100% un¬ 
converted transmission of all orthogonal modes. The 
principal figure of merit for such a drop device is the 
‘selectivity’ [19], defined as: 


S = 


Ipi 


£ N 2 

3 = 1 


(6) 


The best previous attempts at optimizing selectiv¬ 
ity using a single-stage QFC scheme converged on fre¬ 
quency choices that are group-velocity matched, namely 
Pp = Ps 7 ^ P'r( = /?g)[17, 19, 34]. For FWM, this condition 
is easily accessed by situating the channel frequencies on 
either side of the zero-dispersion wavelength of the fiber, 
which automatically grants phase matching for a tun¬ 
able set of frequencies [13]. In addition, for optimum 
selectivity results in FWM, the medium is required to 
be long enough for a complete inter-pump-pulse collision 
(no overlap to no overlap) [34]. However, both TWM and 
FWM schemes have encountered limits that are universal 
to generic systems governed by coupled-mode equations 
such as Eq. (3), and have yielded selectivities limited to 
around 0.8, an example of which is shown in Fig. 2. The 
single-stage selectivity limit arises from in the interaction 
between pulses convecting through each other at differ¬ 
ent speeds. The nonlinear wave-mixing is perfectly phase 
matched only at the central channel frequencies, and the 
phase mismatch acquires opposite relative signs for spec¬ 
tral components of the pulses below and above the central 
frequencies. This causes the QFC contributions of the 
two halves of the spectrum to beat relative to each other, 
which manifest as non-separable time-domain oscillations 
in the GF kernels [19] similar to Burnham-Chiao ringing, 
which occurs when a weak, resonant optical pulse propa¬ 
gates through an atomic vapor [36]. At high conversion, 
the Schmidt-mode profiles become temporally skewed rel¬ 
ative to the group-velocity-matched pump shapes. These 
distortions are minimal at lower CE, enabling our two- 
stage scheme. 


IV. IV. TEMPORAL-MODE 
INTERFEROMETRY 

Our scheme, which we call temporal-mode interfer¬ 
ometry, uses cascaded stages of QFC, as illustrated in 
Fig. 3. The method provides close to 100% TM selectiv¬ 
ity, opening the door for manipulation of TM qubits or 
multi-level qudits, and may also be used for TM chan¬ 
nel multiplexing in classical optical telecommunications. 
In the figure we color code the participating frequency 
channels s in green, and r in blue for the signal pho¬ 
tons, as well as p (and q ) the strong pump field(s) in or¬ 
ange (and red). Our method exploits the fact that using 
single-stage QFC, it is possible to discriminate orthog¬ 
onal TMs nearly perfectly up to a conversion efficiency 




FIG. 2. Single-stage TWM (a,b) and FWM (c,d) dominant 
Schmidt-mode conversion efficiencies \p n \ 2 (a,c) for various 
pump powers, with selectivities S listed in the legend. (b,d) s- 
channel input Schmidt modes, illustrating increased temporal 
skewness with increasing 7 . The values of t' are relative to a 
| P' r — /3' s \l of (b) 20 and (d) 10, where l is the medium length. 
Pump-p is Gaussian in shape. 


(CE) of about 50% [17, 19]. This motivates a two-stage 
interferometric scheme, in which each stage is configured 
for 50% CE for the target TM, and functions as a 50/50 
beam splitter with the r- and 5-frequency channels repre¬ 
senting its two input and output arms. An 5-input pho¬ 
ton in the target TM will be 50% frequency converted in 
the first stage into the r channel with a phase picked up 
from the pump field(s). If all fields are allowed to par¬ 
ticipate in QFC in the second stage with the right rela¬ 
tive phases, the two effective beam splitters will function 
as a frequency-shifting Mach-Zehnder (or Ramsey) in¬ 
terferometer, allowing for complete forward or backward 
conversion of the state of the photon, depending on the 
value of the additional relative phase 0 introduced be¬ 
tween the stages. Between the two stages, the target TM 
component of the signal photon is in a multicolor super¬ 
position of the signal-idler frequency bands (i.e. a color 
qubit) [37, 38]. This interferometric scheme was first pub¬ 
lished by us in a previous work [16]. Phase coherence be¬ 
tween frequency-converted and unconverted components 
of light has since been experimentally demonstrated for 
weak-coherent states by Clemmen et al [26], who inde¬ 
pendently arrived at the frequency-conversion interfer¬ 
ometry idea. However, their implementation used con¬ 
tinuous pumps and so is not TM selective. In our case, 
the interferometric frequency-conversion effect operates 
only on the target TM, as the orthogonal modes have 
negligible CE in both stages. 

There are two configurations of interest: 1) ‘reversed 
collision’ (RC), in which the dispersion in the second- 
stage is inverted relative to that in the first stage, such 
that the relative group velocities of the pulses are re¬ 
versed, and 2) ‘double collision’ (DC), in which the dis- 
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FIG. 3. Temporal-mode interferometer using two nonlinear media (QFC 1 and QFC 2) with identical (DC) or opposite-sign 
(RC) dispersion. Appropriate choices for pump-pulse shapes, pump powers and the phase shift 6 will selectively frequency 
convert a specific s (green, cj s ) TM into an r (blue, uo r ) TM at a different central frequency, while not affecting temporally- 
orthogonal s-input TMs. The pump q is used only for x ^ implementations. WDM stands for wavelength-division multiplexer. 
PC stands for pre-chirp modules, which are necessary for x ^ implementations. The coupler C contains frequency dependent 
delays for the DC case. Reproduced from [16], with permission. 



FIG. 4. Temporal-mode interferometry with the two interferometer ‘arms’ being the frequency channels s (green, lighter shade) 
and r (blue, darker shade), undergoing two complete collisions in the two fibers and exchanging energy. This visualization 
plots the color-coded signal-field intensities in the average velocity frame with populated green input and empty blue input at 
z = 0, for the RC configuration found by numerical solution of the equations of motion, (a) Both pumps are Gaussian, and 
the first Schmidt-mode green-to-blue conversion is nearly 100% at phase-shift 6 = 0. (b) Both pumps are Gaussian, and the 
second Schmidt-mode conversion is nearly 0%. (c) Both pumps are Gaussian, and the first Schmidt-mode conversion efficiency 
is suppressed to zero at phase-shift 6 = n. (d) Pump q is Gaussian, and the shape of pump p is tailored to convert the green 
input mode from (b) into a Gaussian-like blue output with nearly 100% efficiency at 0 = 0. 


persion in the second-stage is identical to that in the first. 
In the DC case, the fast pulses must be time delayed rel¬ 
ative to the slow pulses in between the stages so that 
they walk through each other again in the second stage. 
The inset in Fig. 3 shows the inverse group velocity vs. 
frequency for the two configurations. Figure 4 shows an 
example of TMI for the RC configuration. The neces¬ 
sary time delays for the DC case are implemented in the 
‘coupler’ labeled C in the figure. After the target TM 
is frequency converted, it would be separated from the 
main beam using a standard wavelength-division multi¬ 
plexer (WDM). As shown below, for TWM we predict 
in the RC case a TM selectivity of 0.9846, and in the 
DC case the selectivity is 0.9805, far higher than can be 
achieved with single-stage systems. High TM selectivity 
can also be achieved using FWM, but will require chirp 
pre-compensation of the pump pulses to be implemented 
in modules labeled PC, as described below. 


A. Two-stage theory 

For the two-stage temporal-mode interferometry, the 
combined Green function kernel G rs (t,t f ) is given by the 
interferometric equation [19]: 

/ oo 

&t" [Gr? (t t")G$ (t", t') 

-oo (7) 

T e %e G$ (£, t")G^} (t", t ')], 

where the superscripts indicate the FC process in stages 
1 and 2. Functions and parameters lacking a super¬ 
script stage number characterize the combined two-stage 
process as a whole. For each stage separately, the 
‘transmission’ (no frequency change) coefficient and 
the ‘reflection’ (frequency change) coefficient pn \ both 
taken to be real without loss of generality, indepen¬ 
dently obey relations analogous to beam-splitter rela- 
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tions, Tn ^ + pn = 1- The stage-labeled input modes 
are represented by lower-case functions and 

output modes by upper case <h^). 

For TWM and FWM (with complete pump collisions), 
described by Eq. (3), the Schmidt coefficients are inde¬ 
pendent of pump-pulse shapes, and are determined by 
the value of 7 . In contrast, the Schmidt mode shapes are 
determined by the pump-pulse shapes and the value of 
7 [17]. For FWM, this allows one to change the output 
signal shape to that of any temporal mode, if desired. 
For TWM, which uses of a single pump, one can influ¬ 
ence only the Schmidt mode shapes of the group-velocity- 
matched signal channel. 

The unitarity of the single-stage transformation im¬ 
poses a pairing between the Schmidt modes across the r 
and s channels for that stage [35]. If the input fields for 
a given stage are represented as 

^(OL = ^2 a n^ { nXt'), (8a) 

n 

^(01 in = ( 8b ) 

n 

then the output fields are expressed as 

A P Wlout = + Pn )b n)^nHt), (9a) 

n 

^Wlout = “ Pn )a n)$n\t). (9b) 

n 

The expressions in brackets are equivalent to a beam¬ 
splitter transformation, explaining why the QFC process 
is considered background-free in principle [35]. 

The operating principle of the TMI can be summa¬ 
rized simply as follows. Consider the case that the input 
field to stage 1 is a single temporal mode A^\t) |i n = 
b n <t>n\t)> and the r-input field is empty. Then, from 
Eqs. (4) and (5), the output fields of stage 1, and thus 
the input fields of stage 2 , are 


^Wlout = T n )b n^n ) (t) = A ^ W|i„’ ( 10a ) 

4 1 } (i)| out = P^K^it) = e- ie A<?\t)\. n , ( 10 b) 

where a phase shift 0 of the r field has been introduced 
intentionally by the experimenter before the fields enter 
stage 2. Then the output of stage 2 will be 


A s 2) Wlout =6 « F 

m= 1 

OO 

^Wlout =hn F 

m= 1 


( 11 a) 

(lib) 


where the ‘inter-stage mode overlaps’ are defined as 


Pm,n = J 

( 12 a) 

Vm,n = [ dt^*(t)¥p(t). 

( 12 b) 


If \Pn,n\ = \Vn,n\ = 1, then we say the processes in 
the two stages are temporally mode matched. This oc¬ 
curs only if the dominant output modes of stage 1 co¬ 
incide with the corresponding input modes of stage 2 in 
each frequency channel, that is ^f)n\t ), and 

'J ’n\t) = 

For the target mode n — 1, we wish to have 

A ^(^)lout = 0and7 2 ) (f)| out = & 1 < 7 2 7 ). This can only 
occur if several conditions are met: 1 ) the non-dominant 
Schmidt coefficients are nearly zero; 2) the dom¬ 

inant processes in the two stages are temporally mode 
matched; and 3) the dominant Schmidt coefficients are 
p^ = r-P = y/l/2. Then the phase 0 needs to be ad¬ 
justed to zero (or some other value if the mode overlaps 
are complex). This gives (rpV^/iqi —e ld p^p^rj jp) = 

0 and (pj 2 Vj ^/iqx + e ie r[ 2 ^ p^p ip) = 1. If these condi¬ 
tions are met, then by varying the phase of the r signal 
field between the fibers, either frequency channel in mode 
n = 1 can be 100 % populated at the output of fiber 2 . 

As shown in the Numerical results section, the remain¬ 
ing (non-target) modes (n % 1 ) have p\P ~ 0 , ~ 1 . 

So, light in any one of these modes is not significantly fre¬ 
quency converted in either stage, although there would 
be a change of temporal mode shape if 
creating temporal mode distortion. 

For many applications, especially those involving cas¬ 
caded TM operations, it would be highly beneficial to 
implement TMI without temporal mode distortion of 
the non-converted signals. This is achieved if the out¬ 
put Schmidt modes of stage 1 match the input Schmidt 
modes of stage 2 , that is &n\t) = (j)n\t), and ^n\t) = 
\ipn ; ( t ) for every n. The RC configuration satisfies this re¬ 
quirement along with the conditions §n\t) = <pn\t) = 

$n(t) = 4>n{t ) and &n\t) = 1pn\t) = ^„(i) = 1p n (t) 

( i.e . the output Schmidt modes are identical to the input 
Schmidt modes for both signal channels), leading to the 
elimination of temporal distortion if the non-converted 
input modes are identical to the input Schmidt modes of 
the system. 


B. Numerical results 

In order to derive the Green functions for QFC 
numerically, we implemented a coupled-mode-equation 
solver that accepts arbitrary input functions (A r ( 0 ,t'), 
A s ( 0 ,t')) as arguments, and computes the resultant out¬ 
put functions (A r (/,t), A s (/,t)) for Eq. 3. This is 
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achieved using a fourth-order Runge-Kutta method. The 
solver iterates over differential steps in pulse-propagation 
(A z) from z = 0 to / (medium length). Every iteration 
consists of an upwinded z -propagation scheme for all four 
pulses by a step (A z), followed by a Runge-Kutta im¬ 
plementation of the coupled nonlinear interaction, all in 
spacetime domain. We compute the GF by computing 
the outputs for an orthogonal set of input ‘test signals’. 
To ellaborate, consider the GF kernel and its Schmidt 
decomposition: 

Grs(t,t') = ^p j '$ j (t)4>*(t'). (13) 

3 

The objective is to calculate all the individual compo¬ 
nents (pj, Tj(t) , 0 n (^)) on right-hand side. We first 
pick two arbitrary spanning sets of basis functions { B r 
and {B s g} and re-express the Schmidt modes: 


*j(t) = E u jk BrAt); = E; 

k l 

(14) 


G rs (t, t f ) — 


k,l 


UjkPjVji 




= ^[Grs] kl B r At)Bl l (t')‘ 

k,l 


(15) 


Using (A r (0,£) = 0, A s (0,t') = B s g(t')) as inputs 
for the solver, and decomposing the resulting r-channel 
outputs A r (l,t ) in the {B r , k } basis, will yield the en¬ 
tire Zth column of the complex matrix G rs . Once this 
matrix is determined, its Schmidt decomposition will di¬ 
rectly reveal {Ujk}, {Vji }, and {py}, and through them, 
the Schmidt modes. For all the results presented, we 
chose Hermite-Gaussian functions for the spanning set 
of basis functions for both input and output Schmidt 
modes, since the low-CE Schmidt modes for Gaussian 
pump shapes are nearly Hermite-Gaussian [18, 39]. 


1. Three-wave mixing 

Consider the case in which the target mode is 50% con¬ 
verted in each stage individually while orthogonal modes 
remain almost completely unconverted. Then, using Eqs. 
(5) in Eq. (7) shows the need for inter-stage temporal 
mode matching between the output Schmidt modes of 
the first stage and the input Schmidt modes of the second 
stage for the scheme to work. Figure 5 shows the plots 
for the input and output Schmidt modes of both stages 
for both RC and DC configuration TMI, with a Gaussian 
pump, and the parameter ( = \/3' r — /?' \ l/r p = 200, where 
l is the per stage medium length and r p is the Gaussian 
pump width. The Schmidt modes of a single channel 
for a given stage are found to be temporally symmet¬ 
rically skewed in opposite directions between the input 


and the output, the direction being dependent on the 
sign of (/?' —/?'). Consequently, the inter-stage mode 
matching will be exact for the RC configuration, and in¬ 
ferior for the DC configuration. For ( = 200, the RC 
fjTW M _ rjTWM _ i But f or p>C, = 0.983 and 

rjrwM = 0 . 901 . 

Figures 6(a-d) Schmidt modes for the RC configura¬ 
tion TWM-TMI for ( = 200 and a Gaussian pump. The 
selectivity S was computed to be 0.9846 (\pi\ 2 = 0.9975, 
\p 2 1 2 = 0.0110). (is also the ratio of time-widths of the r- 
and 5-Schmidt modes. The CE of the first four Schmidt 
modes (i.e. \p n \ 2 ) are plotted in green in Fig. 7(a), along 
with the CE of the same Schmidt modes when using a 
custom pump shape (magenta) tailored to ‘drop’ the sec¬ 
ond r-input Schmidt mode. Since the TWM pump can 
influence only the shape of the r-channel Schmidt modes, 
the tailored pump in Fig. 7(b) will frequency convert the 
second r-input Schmidt mode from Fig. 6(a) into an s- 
channel TM identical to the first 5-output Schmidt mode 
from Fig. 6(d). 

As an example, with a Gaussian pump, a signal photon 
in the r channel with shape ijji (t) (Fig. 6(a)) will be 
frequency converted into an 5-channel photon with shape 
4>i (t) (Fig. 6(d)), with an efficiency of |pi| 2 = 0.9975. In 
contrast, a photon in the r channel with shape ^(t) will 
be frequency converted into an 5-channel photon with 

shape 4>2(t), with a very small efficiency of |/? 2 1 2 = 0.0110. 

In other words, an r-channel photon stays in the r channel 
with probability 0.9890, and will exit the device with 
shape T 2 (t) (Fig. 6(b)). Figure 6(c) comes into play 
only if the input photon is in the 5 channel. 

For the DC configuration with ( = 200, we computed a 
selectivity of 0.9805 (with |/q| 2 = 0.9957, |p 2 1 2 = 0.0134). 
This is slightly lower than the RC configuration due to 
relatively inferior inter-stage mode-matching. 

Both single-stage Green function separability, and 
inter-stage mode matching improve asymptotically with 
increasing (, providing for corresponding gains in selec¬ 
tivity (Fig. 8). Larger ( also implies a decrease in tempo¬ 
ral skewness of the Schmidt modes relative to the pump 
shapes. The exact mode matching in RC configuration 
results in matching skewness directions for input and out¬ 
put Schmidt modes of a single channel, a feature that is 
not present for DC configuration (Fig. 9). Therefore, un¬ 
converted higher-order Schmidt modes of the RC process 
undergo no temporal distortion upon passing through 
the TMI transformation. This allows TWM-TMI de¬ 
vices operating in the RC configuration to be used in 
a chained sequence to implement multiple operations on 
the temporal-mode basis. 

TMI is physically easier to implement in the DC config¬ 
uration. A ( of 200 can be realized, for example, in a typ¬ 
ical 5-cm long periodically poled lithium niobate waveg¬ 
uide with a 70-fs pump-pulse and signal wavelengths at 
795 nm and 1324 nm. RC configuration implementations 
rely on media with negative dispersion, which exist for 
y( 3 )-media in the form of photonic-crystal fibers, and are 
also possible to engineer in -media using photonic 
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FIG. 5. The dominant input (a, c, e, g) and output (b, d, f, h) Schmidt modes for the first-(a, b, c, d) and second-(e, f, g, 
h) stages for both r (a, b, e, f) and s (c, d, g, h) channels for both RC (red, lighter shade) and DC (black) configurations of 
TWM-TMI for a Gaussian-shaped pump, and £ = \(3' r — /3' s \l/r p = 200. The values of t are relative to a \/3' r — /3' s \l of 20. Due to 
the nature of temporal skewing, the inter-stage mode-matching between the first-stage output Schmidt modes (b, d) and the 
second-stage input Schmidt modes (e, g) is larger for the RC than the DC configuration, thus yielding better selectivity. The 
complete two-stage composite system Schmidt modes for TWM-TMI in the RC configuration are plotted in Fig. 6. 
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FIG. 6. Schmidt modes for TWM-TMI in the RC configuration with Gaussian pump, and £ = \/3' r — /3' s \1/t p = 200, yielding a 
selectivity of 0.9846. The magnitudes of the first two Schmidt modes are shown for r-input (a), r-output (b), s -input (c), and 
5-output (d). Values of t are relative to a \/3' r — /3' s \l of 20. Fig. 7 shows the conversion efficiencies of the first four Schmidt 
modes for two different pump shapes. 



FIG. 7. (a) Conversion efficiencies for Gaussian-pumped- 

TWM-TMI Schmidt modes in the RC configuration. The 
CE for a Gaussian pump are shown in green (darker shade), 
and those for a custom pump shape tailored to drop the sec¬ 
ond Schmidt mode are shown in magenta (lighter shade), (b) 
The corresponding pump shapes shown with matching colors. 
Values of t are relative to a \/3' r — /3' s \l of 20. 


crystal waveguides [40]. 

TMI can also be extended to arbitrarily large num¬ 
ber of stages in both RC and DC configurations. Since 
the single-stage Green functions are more separable at 
lower conversion efficiencies, increasing the number of 
stages can increase selectivity. Multistage implementa¬ 
tions will use lower pump powers due to the increase 
in the number of inter-pulse interactions. For optimal 
selectivity when using an Wstage process, the conver¬ 
sion efficiency for the first stage should be approximately 
0.5[1 — cos( 7 t/N)]. The expression is the exact splitting 
ratio necessary for a sequence of N nonpolarizing beam 
splitters to interferometrically change the propagation di¬ 
rection of an incident beam of coherent light, and be¬ 
comes more accurate for TMI for larger N. Figure 10 
shows the ratio of converted-to-total energy of the first 
Schmidt mode versus propagation distance for various 
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FIG. 8. Dominant Schmidt mode conversion efficiencies for 
TWM-TMI in the RC configuration, illustrating asymptotic 
improvement in selectivity with decreasing pump-pulse width 
Tp relative to interaction time (/3' r — 


FIG. 10. Multistage TWM implementation of TMI, showing 
the ratio of converted to total energy of the first Schmidt 
mode versus z (propagation distance). N is the number of 
stages and l is the medium length for each stage. 



FIG. 9. The first Schmidt modes of (a) r-input, (b) r-output, 
(c) s- input, and (d) 5-output channels for both DC (black) 
and RC (red, lighter shade) configurations of a two-stage TMI 
implemented using TWM with a single Gaussian pump, and 
\/3' r — /3' s \1/t p = 10. Values of t are relative to a \f3' r — /3' s \l of 
10 . 


number of stages. The plots tend to asymptotically con¬ 
verge to a sinusoidal curve. A four-stage scheme yielded 
DC and RC selectivities of 0.9977 and 0.9978 respectively. 
The corresponding ten-stage selectivities are 0.99996 and 
0.99997. Increasing the number of stages also decreases 
the temporal skewness of the first Schmidt mode relative 
to the corresponding pump-shape. The directionality of 
skewness for 7V-stage TWM-TMI in the DC configuration 
is independent of N. However, for 7V-stage TWM-TMI 
in the RC configuration, one can choose between distor¬ 
tionless FC and distortionless unconverted transmission 
of pulses by choosing odd or even N respectively. 



FIG. 11. The first-stage (a, c, e) and second-stage (b, d, 
f) pump chirps used to derive the three inter-stage mode- 
matched Schmid-mode phase profiles shown in Fig. 12. The 
chirp-parameters {e p , e q } are {1, 0} in (a, b) , {2, 0} in (c, d) 
, and {0.5, 0.5} in (e, f). Values of t are relative to a \/3' r — /3' s \l 
of 10. 


2. Four-wave mixing 

In FWM, there are two strong pulsed pumps, each of 
which needs to be group-velocity-matched to a different 
signal channel to ensure good single-stage Green function 
separability. In addition, the medium being used 
in each stage must be long enough for complete inter- 
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FIG. 12. The dominant input (a, c, e, g) and output (b, d, f, h) Schmidt modes for the first-(a, b, c, d) and second-(e, f, g, h) 
stages for both r (a, b, e, f) and s (c, d, g, h) channels for RC configuration of FWM-TMI for a Gaussian-shaped pumps, and 
complete pump collisions (\/3' r — /3' s \l/(r p + r q ) = 5). Also shown are the corresponding Schmidt mode phase profiles (i-p) for 
three seperate choices of pump chirps for the two stages. Due to pump-induced self- and cross-phase modulation, the second- 
stage pump chirps have to be related to the first-stage pump chirps in the manner specified in Eq. 16 for good inter-stage 
mode-matching between (b, d, j, 1) and (e, g, m, o) respectively. The complete two-stage composite system Schmidt modes for 
FWM-TMI in the RC configuration are plotted in Fig. 13. Values of t are relative to a \/3' r — /3' s \l of 10. 


pump-pulse collision to occur. This is straightforward to 
satisfy with the use of highly-nonlinear optical fibers. For 
Gaussian pumps of temporal widths of the order of 10 ps, 
and wavelengths at 800 nm and 850 nm, a typical highly- 
nonlinear photonic crystal fiber will need to be about 20 
m long for complete collision (| /3' r — /3' s \l/(r p + r q ) = 5). 

The presence of two pumps enables independent shap¬ 
ing of the Schmidt modes of both channels in FWM-TMI. 


However, nonlinear phase modulation (Eq. 3) severely 
affects the Schmidt-mode phase-profiles, and restricts 
FWM-TMI to the RC configuration. This is a result of 
the well-known cross- and self-phase modulation present 
in any medium. To overcome this impairment, the 
pumps must be pre-chirped with specific phase-profiles 
[34] for each stage to enhance inter-stage mode-matching. 
The starting phase profiles of the pump p and q in stage-£ 
( ctp^ ( t ) and ( t ) respectively) need to be 
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(16) 


where 7 ^) = 7 — /3^), and Aj(£) = Aj(0, £). The 
limits of the integration terms are valid only for complete 
pump collisions. It must be noted that the pump phase 
profiles contain a term each that is fully linear in time, 
which are equivalent to frequency shifts. These are an 
approximate compensation for a non-separable part of 
the phase profile of the GF kernels [34] which only attains 
significant magnitude away from the centroids of the GF- 
kernel amplitude functions. The slope of this term for 
each pump is proportional to the magnitude-square of 
the amplitude of the other pump’s envelope. The slope 
is also proportional to 7 ^, which changes sign between 
the two stages. For a given set of chirp parameters, a 
pump will require frequency shifts of opposite signs in the 
two stages. Figure 11 shows the plots of the pump phase 
profiles from Eq. 16 for three sets of chirp parameter 
values. Figure 12 shows the amplitudes and phase profiles 
of the first input and first output Schmidt modes for both 
stages of FWM-TMI in the RC configuration for the same 
three sets of chirp parameter values. 

The chirp parameters e p and e q can be any real val¬ 
ues. However, for FWM-TMI in the RC configuration to 
be useful, one will have to pick specific values to make 
the desired channel input/output Schmidt mode have a 
flat phase profile. For example, Figs. 13(a)-13(d) show 
the Schmidt-mode amplitudes, and phases corresponding 
to FWM-TMI with selectivity 0.9873 (|pi | 2 = 0.9973, 
|p 2 1 2 = 0.0082). The pump pre-chirps (Fig. ll(c, d)) 
were specifically chosen to yield flat phase profiles for the 
r-output (Fig. 13(f)) and 5 -input (Fig. 13(g)) Schmidt 
modes. Namely, e p = 2 and e q = 0. The results were 
computed for |/?' — /3' s \l/(r p + r q ) = 5. 

FWM-TMI can in principle be extended to multiple 
stages, provided that every stage-interface is in the RC 
configuration. The pump frequency shift will alternate in 
sign for every stage in sequence. The pump chirp profile 
functions in Eq. 16 will likewise need to be generalized 
for arbitrary number of stages. The Kronecker -^2 term 
will pick up a factor inversely proportional to the number 
of stages. 

The need for pump-phase chirps in FWM-TMI can be 
circumvented if one resorts to asymmetrically-pumped 
Bragg scattering [41], in which one of the pumps is made 
very long or continuous wave (CW), and the shorter 
pump power is relatively weak. Then, all self- and 
cross-phase modulation effects occur only due to the 
long CW pump, which is equivalent to pump-power- 


dependent frequency and wavenumber shifts of the sig¬ 
nals. This method can simulate TWM-like dynamics in 
a X < 3 ) medium, and by extension, can be used to im¬ 
plement TWM-TMI using FWM. Pump-signal group- 
velocity matching is easier to achieve in media with 
waves whose interaction is phase matched, and the choice 
of which pump is to be made CW does not affect the 
selectivity, making asymmetrically pumped Bragg scat¬ 
tering a valuable option warranting further study. 

V. V. CONCLUDING REMARKS 

TMI can nearly perfectly separate field-orthogonal 
temporal modes into separate beams corresponding to 
different carrier frequencies, as well as combine such TMs 
into a single beam. It also enables real-time routing ca¬ 
pabilities in all-optical networks while transcending the 
limitations of narrow-band filtering and projective mea¬ 
surement schemes [42]. It offers true field-orthogonal re¬ 
alizations of both classical [43] and quantum [42] optical 
code-division multiple-access (OCDMA) networks with 
high spectral efficiency for both channel routing and de¬ 
tection, unhindered by multiuser interference. 

Breaking the 80% single-stage selectivity limit provides 
a means of physically realising high-fidelity quantum¬ 
unitary devices such as the proposed ‘quantum pulse 
gate’ (orthogonal-field TM sorter) and ‘quantum pulse 
shaper’ [44], which operate explicitly on the TM Hilbert- 
space. TMI can aid in long-distance quantum com¬ 
munication and quantum key distribution using novel 
continuous-variable encodings of qubits or qudits in the 
time-frequency basis [45]. The tunable interferometric 
phase-angle allows TMI to simulate linear optical ele¬ 
ments analogous to the action of polarizing and non¬ 
polarizing beam splitters, half-wave plates and quarter- 
wave plates on the polarization basis. TMI devices can 
therefore perform arbitrary unitary transformations on 
the TM space, and in conjunction with multimode quan¬ 
tum memories (such as atomic Raman memories [46, 47]), 
open a new avenue for linear optical quantum computing 
in a single spatial mode [48]. All of these applications 
will be ellaborated upon in future publications. 

In conclusion, TMI can selectively manipulate pho¬ 
tonic temporal mode components of any shape accessible 
by our ability to reshape strong pump pulses [49, 50]. 
Modular TMI devices also possess a phase control (0) 
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FIG. 13. The first two Schmidt modes (a, b, c, d), and first Schmidt-mode phase profiles (e, f, g, h) for FWM-TMI in the 
RC configuration with Gaussian pumps. The first two Schmidt modes are shown for r-input (a), r-output (b), 5 -input (c), and 
5 -output (d). Pump pre-chirp parameters used were e p = 2, e q — 0. The selectivity was 0.9873. Values of t are relative to a 
\P' r -P' 8 \l of 10 . 


that can suppress QFC completely and on demand. TMI 
may be implemented using either three-wave mixing, or 
four-wave mixing enabled frequency conversion, and may 
be generalized to multiple stages for gains in selectiv¬ 
ity asymptotically approaching unity. The addition of 
temporal-mode analysis to polarization and orbital an¬ 
gular momentum analysis completes the toolkit for en¬ 


coding information on quantum and classical states of 
light. 
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